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Abstract 

Motivated by the recent experimental studies on layered ferromagnetic metallic system Gdl2 and its 
doped variant Gdl2Hi; we develop a model to understand their ground state magnetic phase diagram. 
Based on first principle electronic structure calculations we write down a phenomenological model and 
solve it under certain approximations to obtain the ground state energy. In the process we work out 
the phase diagram of the correlated double exchange model on a triangular lattice for the specific band 
structure at hand. 
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I. INTRODUCTION 



Layered magnetic systems with coupled charge and spin degrees of freedom have received 
considerable attention recently owing to their rich phase diagram and intriguing transport prop- 
erties. The layered system Gdl2 and its doped variant Gdl2H2,. , recently studied by Felser et al.— 
and Ryazanov et al.^, show significant magnetoresistance'^ of about 60% at room temperature 
with a ferromagnetic (FM) transition temperature close to 300K— . Unlike the magnetoresistive 
manganites or their bilayer counterparts, there is no significant Jahn- Teller coupling in these 
systems and that enables one to study a purely electronic model, coupled to magnetic degrees 
of freedom, without the involvement of the lattice. In addition, Gdl2 is isostructural to the 
well-known dichalcogenides like 2H-TaS2, 2H-TaSe2 with hexagonal layered structure, but do 
not show any superconducting instability - in fact the resistivity rises at low temperatures. 

It has been pointed out recently- that the large magnetoresistance in Gdl2 is primarily due 
to the freezing out of the spin disorder scattering of the conduction electrons in a magnetic field. 
There is indeed evidence for spin fiuctuation coupled to charge degrees of freedom in the ESR 
experimental. The broad resistivity anomaly at Curie temperature (Tc) shifts towards higher 
temperatures with magnetic field. From susceptibility measurements^ a variety of regions with 
different degree of spin ordering have been proposed. The high-temperature ferromagnetic 
phase gives way to regions of frozen or short-range correlated spins at lower temperatures. At 
about 33% hole doping a sliver of spin glass phase appears at low temperatures (T < 20K) 
changing over to a paramagnetic (PM) phase with doping. 

In Gdl2 the rare earth Gd ions have one electron in the (i-band coupled to the localised Gd 4/ 
electrons through ferromagnetic exchange in the configuration Apbd^. The rare earth ions are 
arranged in a hexagonally coordinated layer structure, each Gd ion has six nearest-neighbour 
Gd ions. Each such layer is separated from the next one by two layers of iodine atoms. This 
results in having the inter-layer coupling among the Gd ions considerably weaker than the intra- 
layer one. LDA band structure results indicate a spin splitting of the conduction banc^ - the 
splitting is nearly complete in the LDA+U results. The crystal field in the trigonal prismatic 
layered dichalcogenide is negligible and the 5d orbitals are quite strongly mixed. Band structure 
calculations (see below) indicate that there is one half-filled d-band that crosses the Fermi level. 
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The on-site Coulomb interaction in the rare earth 5d level is generally not very strong. But a 
single, half-filled narrow d-band crossing the Fermi level and low dimensionality of the system 
make this repulsion quite relevant. The presence of short-range spin fiuctuations and spin 
disordered phases already point to competing magnetic exchange interactions in the system. 
In the half-filled single rf-band, even a small Coulomb repulsion generates AF fiuctuations that 
will compete with the FM interaction mediated by the conduction band via the FM f — d 
exchange (similar to the Zener double exchange (DE) mechanism). It is also useful to note 
that the Gdl2 system never shows full saturation moment, predicted by the LDA calculations. 
The ordered state has a moment of 7.33/iB; less than the full saturation moment of SyU^ out 
of which presumably comes from the 4/ core spinsi. Hence the moment coming from 
the bd electron is considerably less than l/i^ indicating that correlation effects (beyond the 
density functional calculations) and coupling of the spin and charge degrees of freedom are 
very relevant. Also there is almost no literature on the study of a correlated DE model on a 
non-bipartite lattice. On a triangular lattice, it is known that the nearest-neighbour Ising AF 
interactions are frustrated leading to a finite ground state entropy. Indeed, the ground state 
of AF Heisenberg spin model or any correlated electronic model on a triangular lattice is not 
known. In the few studies that exist on frustrated itinerant systems'', the geometric frustration 
of the lattice has often been replaced by a random spin exchange. For a review on these and 
related issues, see, e.g., Ref. 0. 



II. LSDA BAND STRUCTURE 



Spin-polarized band structure calculations were performed for the experimental crystal struc- 
ture of Gdl2^ within local spin-density approximation (LSDA) using the LMTO method^ in the 
atomic sphere approximation with the combined correction term taken into account. Gd 4/ 
electrons were assumed to be completely localized and treated as quasi-core states. All seven 
majority spin 4/ states were occupied while the minority spin ones remain empty which models 
the high-spin state of the Gd 4/ shell. Ferromagnetic arrangement of Gd magnetic moments 
was assumed. 

The calculated band structure shown in Fig. ^ agrees well with the results of previous 
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FIG. 1: (Color online) LSDA band structure and the total DOS of Gdl2. Majority and minority spin 
bands are plotted with solid and dashed lines, respectively. 

calculations^. I 5p states are completely filled and separated from partially occupied Gd 5d 
states by a gap of ~2 eV. The latter are split into majority and minority spin subbands by 
strong exchange interaction with the completely spin-polarized 4/ shell, the splitting being of 
order ~0.9 eV. The total width of 5d derived bands is 5 eV. The most striking feature of the 
band structure is that for each spin direction two lowest Gd 5d bands are split off from the 
others by a gap of 1 eV. 

Because of the trigonal local symmetry of a Gd site d orbitals are split into a one- dimensional 
representation ai (d^^) and two two-dimensional representations ei {d^z, dyz) and 62 {d^y, 
d^2_y2). Band plots show the contribution of different Gd bd orbitals to a Bloch wave function 
represented by the size of a dot (Fig. EI). One can glean that the assumptioni of two lowest 
Gd d bands formed predominantly by a d^^ orbital is not confirmed by our calculation. Except 
for the k space region around the T — A high symmetry direction, the contributions of d^2, 
dxyi and dx2_y2 orbitals to the bands in question as well as to the two bands above the gap are 
comparable. Thus, the gap appears as a result of a rather strong hybridization between the 
above orbitals in the ah plane. Four highest Gd d bands are formed by dxz,yz states which are 
shifted towards higher energies due to the trigonal component of the ligand field. They do not 
hybridize with the former three orbitals. The same conclusion on the nature of the gap was 
drawn by Mattheissi^ based on a tight binding fit to the calculated band structure of a layered 
2H-TaSe2 compound which is isostructural to Gdl2. 
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FIG. 2: (Color online) The expanded view of the Gd 5d bands. The size of filled circles is proportional 
to the squared contribution of 5d orbitals of different symmetry to the Bloch wave function at a given 
k point. The coordinates of the high symmetry points are F (0,0,0), M (-^,0,0), K (-^,^,0), and A 
(0,0,0.136) in 27r/a units; a being the lattice constant. 

III. EFFECTIVE MODEL 

It is clear from the band structure calculations that the band crossing the Fermi level is a 
mixture of primarily three d-orbitals, namely 6.^2, d^y and d^^-yi. It is therefore not possible to 
obtain a one-parameter, tight-binding fit of this band with only nearest neighbour hoppingi^. 
We fit this band using a three parameter (to, tnn and tnnn = -0.514, 0.070 and 0.073 eV) tight 
binding one (Fig. 3a). The band width is about leV and we assume the on-site repulsion 
to be close to or less than this. In any case, we treat this as a parameter in our calculation. 
The corresponding density of states (DOS) is shown in Fig. 3b along with the DOS of the 
one-parameter nearest neighbour DOS used by Hanisch et al^^ in their study of the Hubbard 
model on a triangular lattice. The values of / — exchange have been estimated from the 
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FIG. 3: (a) Tight binding fit (line) for tfie 5d band (dotted line) that crosses the Fermi level in Fig. 
1, (b) DOS corresponding to the tight binding band in (a) with the inset showing the DOS for the 
tight binding band considered in the calculations of Hanisch et. alM- (see text) 

above band structure results assuming a mean-field description for the exchange and a classical 
description for the core spin. Assuming further that the spin splitting is entirely due to f — d 
exchange, the value Jj^S = 0.9eV then serves as an upper limit for the value of the exchange 
coupling {S = 7/2 is the 4/ core spin). 

Based on the above arguments, a possible Hamiltonian for the effective degrees of freedom 
can be written as 

<ij> i <ij><J i,cr 

Here J is the exchange interaction between the localised 4f spins Sj {S = 7/2). The second 
term is the f — d FM exchange coupling (S, s are the core and conduction spin, respectively) 
while the third term represents the hopping of conduction electrons. The last term represents 
the Hubbard interaction. In second-order perturbation theory, J ~ —\tij\^/{JfdS + U) is due to 
the virtual exchange of conduction electrons between sites i and j. It turns out to be antifer- 
romagnetic (AF). Its magnitude is small for the parameters discussed above. Away from half 
filling, the DE mechanism favours an FM configuration of the core spins driven by the kinetic 
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energy (KE) of the electrons, while the Coulomb repulsion U favours an antiferromagnetic ar- 
rangement. This leads to competing ground states tuned by filling, with additional richness 
coming from the intrinsic geometric frustration of an AF state defined on a triangular lattice. 
In the following, we keep J = and study the magnetic phase diagram of the model. We treat 
core spins (S=7/2) semi-classically and use mean- field theory for the Coulomb repulsion U. 

The effective 5d band at the Fermi surface in the undoped system [x = 0) is half-filled i.e., 
the electron density n = 1. Hole (electron) doping, therefore, implies n < 1 (n > 1). Doping of 
holes in the Gdl2 is effected by the insertion of hydrogen^ which localises electrons away from 
the band. This is taken care by suitably moving the Fermi energy. In the limit of Jfd = and 
with a nearest-neighbour tight-binding band, this model represents the single band Hubbard 
model on a triangular lattice studied within mean-field theory by several gronp o^^' ' ^^' ' ^'^' ' ^^' ' ^^^^'^ . 
These studies find a number of magnetically ordered phases, i.e., the FM state, the Neel state 
[ferromagnetic rows coupled antiferromagnetically, the ordering wave vector Q = (27r/v^,0)], 
the AFM state (classical state, angle between spins 120°) and the linear spin-density wave state 
(LSDW), containing zig-zag FM chains [Q = (vr, 0)]. For low U there is a region of no long 
range order (paramagnetic phase) even at half-filling, contrary to the square lattice where the 
band is particle-hole symmetric. 

Disorder in this system plays a significant role in determining the resistivity and short- 
range order, particularly at low temperatures. In fact, the resistivity is too high - at lOOK the 
minimum value is about 20flcm, three orders of magnitude larger than the maximum metallic 
resistivity for this system (about 10 mQcm, assuming a cylindrical FS and lattice constants 
a,b c::^ AA along the plane). It cannot be only disorder that is responsible for this high resistivity 
in the metallic state, the interaction surely plays a major role. At the moment we do not take 
the disorder into account and look at the clean limit to investigate the possible magnetic orders 
that the model provides. The effects of disorder will be taken up in a separate calculation^^ 
later. 

It is also useful to note that the 4/ levels are not very far from the Fermi level and a small 
hybridisation Vfd between 4/ and 5d electrons cannot be ruled out completely on the basis of 
energetics alone from the band structure results. However, such a term should be very small 
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FIG. 4: (Color online) Phase diagram in the U — n plane at Jjd = 0, 0.005, 0.07 

owing to reasonably large Uff and at most add a small AF coupling between the d and / 
electrons effectively reducing Jf^ only slightly. 

IV. CALCULATION AND RESULTS 

The physics of the model is controlled by competition between the FM double exchange and 
the AF correlations coming from the Coulomb repulsion U. The geometric frustration of the 
AF state on a triangular lattice adds to the number of possible ground states and their closeness 
in energy. In the limit Jfd —>■ oo, one can project out the direction of conduction electron spins 
parallel to the local core spin by using the transformation e*'^z<^i/2g«o-yf!/2g-«o-z<^i/2_ operates on 
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a two-component spinor and reduces to a system of spinless fermions (the angles 9i , refer to 
the polar and azimuthal angle of the conduction electron spin with respect to Sj). The on-site 
Coulomb term, in that case, becomes irrelevant. In this situation the AF superexchange term 
becomes important, though the DE mechanism would lead to an FM state except close to n = 
and n = 2. 

However, the Gdl2 and Gdl2H2^^ systems are not really described by this limit and a finite 
Jfd is the relevant limit here. We treat the core spins classically and the Coulomb interaction 
term in a mean-field approximation. To describe the different magnetic states in this approxi- 
mation, we choose S = SoC*^ ""', where Q is the wave vector corresponding to the chosen state 
described earlier. We choose the spin quantization axis in the xy plane. The second term in 
the Hamiltonian Eqn.(l) then connects every k point in the Brillouin zone to k ± Q and one 
has to diagonalise the Hamiltonian on a k-grid in the hexagonal BZ. The order parameter is 
obtained by minimising the free energy for a particular filling. For the present calculations 
we have chosen to investigate the following ordered phases known from previous investigations 
on triangular lattices. The FM phase, the Neel state, the three-sublattice AFM state and 
the LSDW state mentioned earlier. We use the three parameter tight-binding fit to the band 
appropriate for Gdl2 and Gdl2Ha; systems for our calculations. 

As mentioned above, the Coulomb interactions are generally not very strong in the bd level, 
though for the narrow band crossing the Fermi level (see Fig.l) they may still be very important. 
Since both f — d exchange and Udd contribute to the spin splitting of the conduction band, it 
is difficult to specify these two parameters. We therefore allow them to vary over a possible 
range, treat them as free parameters and draw the magnetic phase diagram. 

In Fig. 4 (a)-(c) the phase diagrams for Jfd = 0, 0.005 and 0.07eV are shown. The Jfd = 
phase diagram can be compared with a similar phase diagram obtained by Hanisch et al.^^ with 
a nn tight-binding band dispersion = 4t cos ^^^^|^ cos ^ + 2t cos ky. It is noticed that there 
are considerable differences between the two calculations. At half- filling, n = 1, we observe the 
Neel state for U close to OAeV (about 0.36Vr, where W is the bandwidth). This continues all 
the way upto U = leV. The region of Neel state is quite wide in our phase diagram around 
n = 1. For U < 0.39el^, we obtain a paramagnetic phase all the way down. Hanisch et al. 
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report three-sublattice AF phase at half-filling at Ured = u^uc ~ ^'^^ = 15.81t for the 
triangular lattice, so that U = 0.58W^) similar to what has been observed earliei»i^. Away from 
half-filling we get a small region of phase segregation (PS) followed by LSDW, AFM, FM and 
finally another region of PS depending on the value of U on the hole-doped side, while for n > 1 
there is a large region of phase separation from the paramagnetic phase. 

The phase diagram changes dramatically with finite values of Jfd- As seen in Fig. 4(b) the 
PM phase reduces to a region close to the empty (n = 0) and full (n = 2) bands. Due to 
the asymmetry in the DOS, the regions of stability of the PM phase for n = and n = 2 are 
very different, the latter being larger. The rest of the PM region becomes FM for finite f — d 
exchange as expected. The AFM region completely disappears while the Neel region becomes 
narrow beyond U = leV. The LSDW phase also shrinks considerably. With increasing Jf^ 
the FM region takes over completely (Fig. 4(c)). Two narrow PM regions survive at the ends 
owing to negligible KE contributions. 

It is obvious from the phase diagrams shown that the value of Jfd for Gdl2Ha; systems should 
be very small but finite. In the model without disorder, the phases are extremely sensitive to 
changes in the f — d exchange. This is suggestive of the fact that indeed in Gdl2H^ several 
phases appear very close by in energy. Just as in other similar correlated systems with competing 
interactions, e.g., the manganites^^ and the double perovskites^ii, a fine tuning of the model 
parameters is necessary to specify the phases of the system. Such a situation leads to multiple 
phases with possible microscopic phase segregation as we observe in Fig. 4. The appearance of 
a multitude of phase segregated regions in the present case also has to do with the underlying 
frustration of the AF states on a triangular lattice. The large entropy present in the ground 
state or very close to it make the ground states easily tunable. 

In view of such sensitivity of the phases to Jfd, we draw the phase diagram also in n — Jfd 
plane for a fixed value of f/ = O.SeV. Fig. 5 reveals that in order to obtain multiple phases, 
Jfd should be less than about 17.5meV. There are the phases LSDW, Neel, FM and PM 
present in the phase diagram with a region of phase separation coming from the first order 
transition between the Neel and FM phases. For bigger values of Jfd the entire phase diagram 
is ferromagnetic. Such a low value for f — d exchange indicates that the spin splitting in the 
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FIG. 5: (Color online) Phase diagram in the Jjd — n plane for U = 0.8 

d— band is not due entirely to the f — d exchange. Instead a considerable part must come from 
correlations in the conduction band. 

The phase diagram presented in Ryazanov et al.,^ does not seem to show regions of phase 
segregation or AF spin order. On the other hand they do indicate the presence of a spin 
glass region (or a mixture of FM and SG) beyond n = 0.33 as can be inferred from the broad 
susceptibility peak they observe. Doping beyond n ^ 0.7 suppresses long range magnetic order 
completely and leads to a broad paramagnetic region. 

The region ri < 0.7 showed large thermal hysteresis in the susceptibility data. There is 
a splitting between the field-cooled and zero field-cooled susceptibility at low temperature 
indicating disordered magnetic behaviour (arising from domain rotations, magnetic clusters or 
superparamagnetic regions). The ferromagnetic transition temperature is also very sensitive 
to the method of sample preparation and the measurement techniques. As reported by the 
authora^i^ the data are quite noisy and the system appears to be intrinsically disordered, no 
matter how carefully prepared. The high value of resistivity in the metallic phase and the 
unsaturated magnetic moment deep inside the FM phase are suggestive of the presence of 
disorder as well as correlation. All this indicate a possible competition between different ground 
states, high degree of magnetic disorder and a possible microscopic phase segregations. As in 
the observed phase diagram, there is a predominance of the FM region close to half-filling that 
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we also find, but there are other phases close by depending on the value oi f — d exchange 
coupling in the calculated phase diagram. The actual value of hole doping in Gdl2H2; for a 
given value of x, is difficult to ascertain^, though increasing x is understood to have the effect 
of hole doping. It is expected from our calculations that there would be several different phases 
and microscopic phase segregations in this system too. 

The large magnetoresistance in Gdl2 has already been attributed^ to the quenching of spin 
disorder scattering of the conduction electrons. In a microscopically phase segregated system, 
there is likely to be very large magnetic scattering across microscale ordered regions. In an 
aligning field such scattering is reduced drastically leading to a drop in resistivity. In the model, 
we had a superexchange interaction to begin with which we have neglected in the calculations 
that followed. The presence of such a term would increase AF tendencies and thereby allow 
for a larger value of Jfd than we derived above. Besides, the end regions of the phase diagram, 
= 0, 2 will also become antiferromagnetic owing to the absence of itinerant electrons there. 

In conclusion, we have studied a model which, we believe, describes some of the features of 
the nearly two dimensional, triangular lattice magnetic system Gdl2Ha. . The resulting phase 
diagram has some broad similarity with the experimental one, though it predicts much more 
structure in the ground state phases than has so far been observed. We hope to motivate more 
experiments in these directions in order to check on our predictions. 
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